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^ ; ABSTRACT 

Qlh! We test and present the application of the full rescaling method by Angulo & White 

(2010) to change the cosmology of halo catalogues in numerical simulations for cos- 
\^ ' mological parameter search using semi-analytic galaxy properties. We show that a 

^ , reduced form of the method can be applied in small simulations with box side of 

^ ■ ~ 50 h~-'^Mpc or smaller without loss of accuracy. We perform statistical tests on the 

accuracy of the properties of rescaled individual haloes, and also on the rescaled pop- 
ulation as a whole. We find that individual positions and velocities are recovered with 
almost no detectable biases, but with a scatter that increases slightly with the size of 
the simulation box when using the full method. The dispersion in the recovered halo 
[Ci, ■ mass does not seem to depend on the resolution of the simulation. Regardless of the 

' halo mass, the individual accretion histories, spin parameter evolution and fraction of 

l/^ , mass in substructures are remarkably well recovered. In particular, in order to obtain a 

more accurate estimate of the halo virial mass, it was necessary to apply an additional 
I correction due to the change of the virial overdensity and the estimate of its effect on a 

NFW virial mass. The mass of rescaled haloes can be underestimated (overestimated) 
for negative (positive) variations of either as or Qm , in a way that does not depend on 
the halo mass. Statistics of abundances and correlation functions of haloes show also 
small biases of < 10 percent when moving away from the base simulation by up to 2 
times the uncertainty in the WMAP7 cosmological parameters. The merger tree prop- 
^ , erties related to the final galaxy population in haloes also show small biases; the time 

' since the last major merger, the assembly time-scale, and a time-scale related to the 

stellar ages show correlated biases which indicate that the spectral shapes of galaxies 
would only be affected by global age changes of ^ ISOMyr, i.e. relatively small shifts 
in their broad-band colours. We show some of these biases for different separations in 
the cosmological parameters with respect to the desired cosmology so that these can 
be used to estimate the expected accuracy of the resulting halo population. We also 
present a way to construct grids of simulations to provide stable accuracy across the 
Qm VS. as parameter space. 
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1 INTRODUCTION 

The Lambda Cold Dark Matter (ACDM) cosmological 
model is the standard theoretical framework for structure 



formation in the Universe. In order to understand how galax- 
ies form and evolve in this cosmological context, we must 
also understand the properties of dark matter haloes over a 
wide range of physical scales and across the cosmic history. 
Numerical simulations provide one of the best methods for 
approaching this problem. 
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The cosmological parameters that provide the best 
match between the ACDM cosmology and several observa- 
tions, are most often obtained using measurements of the 
power spectrum of tem perature fluctuation s in the cosmic 
microwave background ll Sanchez et al. 20061: Hinshaw el alj 



I2OO9I : iDunklev et al ] I2OO9I : IJarosik el al.l I201W or of den- 
sity fluctimti£ns_^^ distribution I Sanchez et al.l 
I2OO6I : iPercival et al.|[2007h : but these are limited to the lin- 
ear regime of density fluctuations. In order to make com- 
parisons between the model and observations in the non- 
linear regime, without making simpli fled assumptions such 
as th e spherical or ellipsoidal collapse (|Sheth. Mo fc TormenI 
one would need to use fully non-linear numerical sim- 
ulations. However, these are too expensive in terms of com- 
putational time. Great efforts go into running even single 
simulations corresponding to one set of cosmological param- 
eters, which can still be compared with a very wide range of 
observational measurements. 

F or the past several years the Millennium Simula- 
tion (|Springel et al.1 120051 ) has been the focus for many 
studies of the distribution and statistical properties of 
dark matter haloes and provides the basis for the imple- 
mentation o f semi-analytic models of the evolving galaxy 
population (|Bower et al.l I2OO6I : [Croton et al.1 12OO6I '). How- 
ever, it is important to appreciate that this simulation 
was run using the WMAPl cosmological parameter set 
l|Spergel et ahlboOSl ). which are rather different from the 
current best fit parameters. More recently the Bo lshoi Sim- 
ulation (jKlvpin. Truiillo-Gomez fc Primackll2010l ) used cos- 
mological parameters consistent with the l atest WMAP5 
llHinshaw el al.1 l2009l : iKomatsu e t al." 200i; iDunklev et"al] 
l2009h and WMAP7 (jjarosik el al., .2010) values. The main 
difference is that the Millennium simulation used a substan- 
tially larger amplitude of perturbations than the Bolshoi 
one. 

The differences in the parameters of the cosmological 
model can affect the abundances and properties of the dark 
matter haloes at a given redshift. Consequently, the prop- 
erties of the galaxies hosted by such haloes may also be 
affected. As a consequence the interpretation of many obser- 
vational statistical tests, like those measuring clustering and 
luminosity functions, becomes more difficult. Notice that 
according to the cosmological para meter constraints from 
WMAP5 and cluster abundances (|Rozo et al.1 12009| ). the 
Millennium Simulation i s about 2 — a away from the best fit 
model (see Figure 1 in iKlvpin. Truiillo-Gomez fc Primeic^ 
I2OIOI ). The small variations of the parameters of the back- 
ground cosmological model also become important in the 
understanding of the baryonic processes i nvolved in galaxy 
formation. An extreme case can be seen in I Cole et al .1(11994 ) 
where the maximum of the stellar formation rate (SFR) ac- 
tivity is found a,t z — 1 due to their choice of $7 = 1, which 
favours mergers at low redshifts. 

Many of the differences between the Millennium and 
Bolshoi simulations can be re-scaled using different methods. 
A first attem pt to change the co smology of a simulation was 
presented by[Z heng ct al. ( 2002) which was later applied by 
iHarker. Cole fc Jenkins (2007) with the aim of constraining 
the cosmological parameter related to the amplitude of lin- 
ear fluctuations in spheres of 8fe~^Mpc, as, u sing the GAL- 
FORM semi-analytic model (jBaueh et al.ll2005l l and only two 
individual dark matter simulations. This approach already 



greatly diminished the computational time that would have 
involved running numerical simulations for each set of cos- 
mo logical parameters . More re cently, the method suggested 
by lAngulo fc White! (|2010. AWIO ) allows a more flexible 
change in the cosmological parameters of a simulation with- 
out incurring in a important loss of precision. This algorithm 
scales the output of a cosmological TV-body simulation car- 
ried out for one speciflc set of cosmological parameters so 
that it represents the growth of structure in a different cos- 
mology. The accuracy of the rescalings can only be estimated 
at first order using extended Press-Schechter theory or, more 
accurately, by running numerical simulations and testing for 
particular characteristics. 

However, the algorithm developed by AWIO is applied 
to every particle in the simulation and the post-processing of 
the outputs of the simulation must be repeated (i.e. identi- 
fication of the haloes and construction of the merger trees). 
This can be computationally demanding, particularly in the 
case when the rescaling needs to be done several times. We 
are interested in studying whether applying the method to 
dark-matter haloes instead of to particles produces similar 
accuracies in the rescaled simulations, since this would al- 
ready considerably reduce the computational time. 

We study whether the use of the AWIO method allows 
to explore variations in the phenomenology of non-linear 
density fluctuations in the cosmological parameter space 
via, for example, monte-carlo sampling. This is the reason 
why we will concentrate on reducing as much as possible 
the computational time in obtaining a catalogue of haloes 
from an A'^-body simulation on a different cosmology using 
AWIO, and also on obtaining even a reduced version of their 
method. In any of these cases, the computational time in- 
volved is dramatically reduced in comparison to what would 
be needed to run a complete simulation for each new set of 
cosmological parameters. In the case of the reduced algo- 
rithm, it can be applied to small simulations where linear 
corrections span scales of the order of the box size. In the 
full and reduced cases, the approach we follow can be used 
when the particle information is not available as will be the 
case of future large simulations. We perform a number of 
tests using A^-body simulations in similar cosmologies than 
those corresponding to the Millennium and Bolshoi simu- 
lations, in order to quantify the accuracy of the recovery 
of individual and statistical properties of the dark matter 
haloes and also exploring the effect of different extrapola- 
tion baselines. 

The outline of this paper is the following. In Section 
[2] we present a brief description of the rescaling technique 
presented by AWIO. We also describe in this section the 
Ai'-body simulations used in this work, and the aplicabil- 
ity of a reduced version of the method. The results of the 
statistical and individual properties of rescaled dark matter 
haloes are shown in Section [S] In Section |4] we measure how 
the resulting halo catalogues are affected by the distance in 
the cosmological parameter plane Qm — f"8 around the base 
WMAP7-ACDM model; this can later be used to evaluate 
how precise is any rescaled dark matter halo catalogue to 
ensure optimal explorations of the cosmological parameter 
space. Finally, the conclusions are presented in Section [5] 
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Table 1. Relevant parameters used in the simulations. The columns show the name of the simulation, the number of particles, the 
gravitational softening, the initial redshift, the particle mass, the boxsize, the matter density parameter, the baryon density parameter, 
the primordial spectral index, the dimensionless Hubble parameter, and the linear fluctuations amplitude in spheres of 8^~^Mpc. The 
cosmological constant density parameter is Q\ = 1 — Qrn (flat models) in all cases. The boxsize of the B simulations were calculated 
using the resulting scaling factor s obtained from the rescaling technique. 
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2 SCALING THE HALO CATALOGUES OF 
iV-BODY SIMULATIONS 

We now briefly describe the procedure used to scale halo 
catalogues from a given cosmology into a new set of cosmo- 
logical parameters. This procedure consists on the method 
presented by AWIO which can either be applied to haloes or 
individual particles in the simulation. If P{k) is the linear 
matter power spectrum at 2; = we can define the variance 
of the linear density field as, 

\2 



D{z)_ 

4-K 



k^P{k)W\kR)dk, 



(1) 



where ii is a comoving smoothing scale, D{z) is the linear 
growth factor normalised so that D{z = Q) = 1 and W{x) 
is the Fourier transform of a spherical top-hat filter defined 

by, 

sin [x) — X cos {x) 



W{x) 



(2) 



Assuming that we want a halo catalogue of a given cosmol- 
ogy (denoted with B) evolved to a final redshift Zg starting 
from another halo catalogue which has different cosmologi- 
cal parameters (denoted with A), the procedure is to find a 
length scaling s of the boxsize, and a final redshift z;^ defined 
so that the linear fluctuation amplitude ^^(s^^T?, z;^) over 
the range [s~^ Ri , R2\ is as close as posible to aB{R,Zg) 
over the range [Ri,R2]. This is done by minimising the func- 
tion 
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R ' 



(3) 



over s and z'^. The values of Ri and R2 are selected so that 
M{s~^ R2) is the mass of the largest halo in the original 
simulation and M{s~^Ri) the mass of the smallest one. 

Once we have the two parameters s and z^, the boxsize 
is scaled so that Lb = sLa and the earlier redshifts in the 



new cosmology {zb < Zg) arc obtained from those in the 
original cosmology {za < z^) through 



Db{zb) = Da{za 



DA{zi)' 



(4) 



After scaling the positions with s, and having matched the 
cosmological times, we must consider that the velocity and 
mass of the haloes need to be corrected using 



{1 + za) Db{zb) hA ^ 
= *'7-;— — : — 71 — '''a, 



(1 + zb) Da{za) hB 



Me 



s'^Ma, 
Ha 



(5) 

(6) 



where the dot indicates a derivate with respect to time, h 
is the dimensionless Hubble parameter and Q.b and CIa are 
the values of the matter density parameter Q,ra ~ + ^dm 
for different cosmologies, being f2(, and fldm the baryon and 
dark matter density parameters, respectively. 

The steps presented in equations (1) to (6) will be re- 
ferred to as the reduced AWIO method. The full implementa- 
tion requires to correct the positions and velocities for resid- 
ual differences in the power spectrum of the two cosmolo- 
gies and/or for the length scaling. To do this it is necessary 
to modify the contribution of the long wavelength compo- 
nents on the position and velocity flelds using the Zel'dovich 
approximation. The range of modes were the correction is 
applied is fe < k„i, where kni satisfies A^(fe„j) = 1 being 
A^(fc) = fc^P(fc)/27r^ (See AWIO for the full details on this 
correction). 



2.1 iV-body simulations 

In order to test the scaling of haloes and their histories we 
use two sets of simulations. In a first set we use 6 simular- 
tions. The two main simulations are simulation A which has 
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Vinal Mass [M^ h ^] Concentration parameter Cg (v^^j/v^^ Concentration parameter Cg (Bullocl^ Concentration parameter Cg (Zliao) 

Figure 2. Comparison of individual properties between rescaled and reference haloes. The top boxes show the relative difference in 
number of particles (A), the difference between halo positions as a function of the number of particles per halo (B), the relative difference 
between the modulus of the velocity vectors as a function of the peculiar velocity of haloes (C) and the histogram of the angle subtended 
by the velocity vectors of matched haloes (D). The bottom boxes sh ow the virial mass (E) and three different halo concen t ration 
parameters, obtained from the ratio Vmax /v^ir (F), computed using the iBuUock et ah I ll200lli recipe (G) and usig the lZhao et al. 1 1|2008| ') 
prescription (H). In all panels we show the full population of haloes. The black solid lines show the median, and black ashed lines enclose 
80 percent of the haloes. We also show the recovered number of particles for the simulation with a large box and with a lower resolution 
(panel A, green and red lines respectively), the median velocity offset for the large simulation (panel C, green solid line) and the median 
of the relative error in the virial masses before the correction of &^ir{z, ^m) via a NFW profile (panel E, dotted line.) 



Table 2. Rescaling parameters. In the first column we show the 
simulation name, the second and third contain the scaling factor s 
and the final redshift z^, respectively, obtained after minimising 
Eq. [3] The last column shows the minimum (S^ms resulting from 
minimising the difference between the rescaled and desired linear 
fluctuation amplitudes. With the exception of simulation Bbig 
where the range [Mi, M2] = [1 X 10^1, 5 X IQl^Jh-iM© is used, in 
all simulations we use the mass range [Mi, M2] = [1 X 10^", 5 X 
10^'^]h~^}AQ to perform the minimisation. 
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a Millennium-like cosmology and simulation B which has 
the background cosmology of the Bolshoi simulation. Both 
simulations contain 256"^ particles in cubic volumes of 60 



and 67.68h~^Mpc of side length, respectively, which results 
in a particle mass comparable to that of the Millennium 
Simulation, Mp ~ lO^h~^M0. This set of simulations also 
includes a low resolution version of simulations A and B, 
with 128^ particles, and versions with volumes ~ 100 times 
bigger (with respect to A and B), which we will refer to as 
Alow, Blow and Abig, Bbig, respectively. With the low res- 
olution simulations we will minimise the rescaling function 
using the same parameters as for the A and B simulations; 
in this way we will only test the effect from the coarser res- 
olution. In the case of the larger volume, the minimisation 
will be performed on a range of masses one order of magni- 
tude higher than for A and B and will therefore include the 
effects from considering different scales in the power spectra, 
as well as a lower resolution. 

Unless otherwise stated, the B cosmology will be taken 
as the reference model and the simulations with the A cos- 
mology will be rescaled to resemble the former as closely as 
possible. 

The second set contains 9 low resolution simulations 
(128'^ particles) designed to measure the variation in the 
precision of the rescaling of the halo catalogues as the de- 
sired cosmology moves further away in the plane flm — os 
from the original parameters. In order to do that, we fix all 
parameters and vary only Qm (or erg) in ila and ±2cr of 
its original value using a ~ 0.03, which is in concordance 
with the standard WMAP7 deviation for both parameters 
(jjarosik el al.ll2010l ). giving us 8 simulations. The remaining 
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Figure 1. Top panel: ID rms difference in the rescaled particle 
positions for the full (solid) and reduced (dashed) AWIO method. 
The lines show the variation in the errors for fixed redshifts, as a 
function of the simulation box size; the pairs of solid and dashed 
lines correspond to different redshifts. Bottom panel: ratio be- 
tween the ID rms difference in the particle positions obtained af- 
ter applying the reduced and full AWIO method; each line shows 
the results when fixing the box size but allowing the redshift to 
vary (box sizes are shown in the figure key). The x-axis shows the 
product of the box side of simulations Lboa; and the non-linear 
limit mode fc„;(z). 



simulation was run varying both thefi™ and as parameters 
by +la. The relevant parameters for both sets of simulations 
are detailed in Table [T] The range of masses over which the 
minimisation is performed is the same as for the A and B 
simulations, in order to mimic the accuracy of this proce- 
dure for a simulation with the particle mass resolution of 
the Millennium Simulation. 

All the simulations were evolved from their intial red- 
shifts using the public version of GADGET2 (iSpringel.2005iy 
The initial conditions we r e con structed using the public code 
GRAFIC2 (jBertschingeil 1200 ih and use exactly the same 
random seed in all cases. The halo catalogues, including 
substructure identification and merger histories, were con- 
st ructed using the SUB FIND algorithm explained in detail 
m ISpringelet aLll|200ll ). The linking length parameter used 
is equal to 0.17 times the mean interparticle separation and 
considering only groups with at least 10 particles. The simu- 
lations outputs consist of 100 steps equally spaced in log(a) 
between z — 20 to z = 0. 

In the process of obtaining the B simulations we search 
for the scaling factor s and final redshift z^ that provide 
the best fit between the actual linear fiuctuation amplitudes 
via the minimisation of Eq (|3]). These are shown in Table 
[2] for each simulation (see the Table caption for the mass 
ranges used in the minimisation). The table also shows the 
minimum value of (J^ma obtained. 



2.2 When is it acceptable to use the reduced 
AWIO method? 

In this subsection we perform a test that will allow us to 
infer in which cases it is acceptable to ignore the quasi- 
linear correction of long wavelength contributions. This test 
consists on implementing the reduced and full algorithm and 
studying the resulting differences as a function of the boxsize 
of the simulation, and of the redshift. The latter is done since 
kni{z) depends on the redshift. 

We use two sets of four 256^ particle simulations specific 
for this subsection with reference simulation boxsizes Ltox ~ 
50, 100, 500 and 1000 h'^Mpc. The first set has the same 
cosmological parameters than simulation A and the second 
one has the parameter set of simulation B. We will analyse 
the z =0, 0.5, 1, 2 and 4 outputs. 

The top panel of Figure [T] shows the ID rms difference in 
the particle positions, for the case where the AWIO method 
is applied in full and in its reduced version (solid and dashed 
lines, respectively); each line shows the change in the error 
as the simulation box increases to the right. As can be seen, 
for simulations of ~ 50/i~^Mpc a side, the reduced rescaling 
algorithm produces positions that are as precise as those 
obtained using the full method. 

The bottom panel shows the ratio between the ID rms 
differences obtained from the reduced and full AWIO meth- 
ods for different values of the simulation box, as the redshift 
increases to the right. In this case it can be clearly seen that 
the small simulation box (with side ~ 2-n /kni{z — 0)) shows 
almost a unit ratio for all the redshifts explored. 

The fact that the large-scale correction is not important 
in simulations of ~ 50/1" ^ Mpc a side is expected since, by 
construction, the displacement fields are smooth on modes 
larger than k„i. Note that this correction depends on the 
cosmologies selected for the simulations, and therefore the 
results shown in this subsection are only presented as a qual- 
itative example. 



3 RESULTS 

In this section we apply the rescaling to the simulations pre- 
sented above, and perform tests on the recovery of the prop- 
erties of individual haloes, including their detailed growth 
histories. We also study the recovery of statistical proper- 
ties of the global population of haloes. From this point on 
we apply the reduced version of the AWIO algorithm to dark 
matter haloes in all the simulations except Bbig, to which 
the full method is applied since has a boxsize for which the 
quasi-linear correction is important. 

3.1 Comparison of Individual Halo Properties 

We compare the properties of haloes of at least 50 particles 
in the numerical simulations run with the desired cosmo- 
logical parameters (B simulations, the reference model) to 
the haloes from simulations A (with a different cosmology) 
rescaled to the cosmology of the B simulations (rescaled 
haloes). These haloes are allowed any number of member 
particles (at least 10). In order to do this we match haloes 
that share the largest percentage of particles via their GAD- 
GET2 identifying particle ID (referred to as matched haloes) ; 
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Figure 3. Comparisons between individual haloes belonging to the rescaled and reference catalogues. The number of particles per halo 
is shown in the upper panels. From top to bottom the panels show the cumulative accretion history, cumulative accretion via mergers, 
cumulative accretion of individual particles, the evolution of the spin parameter and the fraction of mass in substructures, all as a function 
of the expansion parameter a. The solid lines show the true evolution; the dashed lines are for the evolution in the rescaled catalogue. 
The haloes correspond to the 1st, 10th, 100th, and 1000th largest haloes in simulation B and their FoF masses are 4.2 X 10^^, 8.2 X 10^'^, 
1.2 X 10^3 and 1.1 X 10^^ H'^Mq respectively. 



this can be done as the initial conditions are constructed us- 
ing the same random seed. The percentage of haloes in the 
B simulations that have a matched halo in the A simulation 
is higher than 99% in all cases. 

We compare the properties of matched haloes in Figure 
[21 The number of particles of the recovered haloes shows a 
slight underestimation of < 5 percent, and a clear increase 
in the dispersion that decreases with halo mass and becomes 
10 percent for haloes of ~ 200 particles. Results for a lower 
resolution simulation (Blow) in green, and for the simula- 
tion with ~ 100 times larger volume (Bbig) in red are also 
shown. As can be seen, neither the lower resolution nor the 
larger amplitude over which the minimisation is done affect 
the precision of the recovered number of particles per halo 
for haloes with > 200 particles. The catalogue of simulation 
Bbig shows a slightly broader range of recovered number of 
particles for haloes with less particles, but it is not clear 
whether this is a result of noise or the larger box of the sim- 
ulation. The results for the Blow and Bbig simulations are 
also similar to those of the B simulation in most of the com- 
parisons shown in the other panels of this figure. Therefore, 
in order to improve clarity, in the remaining panels we will 
only show the results for the B simulation except for the 



cases where there are noticeable differences with resolution 
or box size. 

As can be seen the positions of haloes are well recovered 
to a precision of 0.1h~^Mpc for the catalogue of simulation 
A (difference between percentiles 10 and 90 and the me- 
dian). The velocities show relative differences of less than 
~ 5 percent, and small changes in the direction of the ve- 
locity vector with a mode of 1 degree. As can be seen in the 
figure, the velocities tend to be biased high for v < 350km/s, 
and biased low for v > 700km/s, although always below a 5 
percent difference. The precision in the recovered halo po- 
sitions and velocities are consistent with those reported by 
AWIO in their Figure 8 using the full implementation of the 
algorithm. 

The larger box (blue line) produces an increase in the 
bias at low peculiar velocities but removes it at the large 
velocity end. Also, the amplitude of the differences between 
velocities increases more rapidly for low velocities in the 
catalogue of the large volume simulation. 

We compute the virial mass of haloes as the mass in- 
side a sphere that contains an average virial overdensity 
of 5„i,.(z,n„) ^ 1871^(1 + 0.4093a:^'^^'^^ )(l + x^)~\ where 
X = (l/f7m(z = 0) - 1)1/3(1 -I- z)-i (|Nakamura fc Sutol 
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1 19971 ). Since the methodology we apply does not involve re- 
identifying the haloes using dark-matter particles, the virial 
mass for the rescaled haloes resulting from applying Equa- 
tion ([6]) is underestimated by a ~ 10 percent (panel E, dot- 
ted line). However, by considering a correction due to the 
difference in 5„ir(z,f2m) which is 15 percent higher for the 
cosmology of the recovered haloes, we are able to reduce 
the discrepancy in the virial masses to below 5 percent. To 
do this correction we assume a NFW dark matter density 
profile l|Navarro. Frenk fc Whit"3ll997l ') for each halo, com- 
pute the integral of the assumed density profile above the 
corresponding virial overdensity i n each cosmology us ing a 
concentration parameter given bv lBuUock et all (|200lh , and 
then multiply the rescaled mass by the ratio between the 
integrated mass in the target and base models. 

We also explored the effect of adopting different con- 
centrations on the recovered virial masses (Figure [2|. We 
used three different concentrations, (i) the proxy Vmax/vvir 
(where Vmax is the maximum circular velocity of the halo 
and v^ir is the circular velocit y at th e virial radius), (ii) the 
recipe given by 'Bullock et alj (|200ll ). and (iii) from the fit 
provided by i Zhao et al. (2008h . Our results indicate a negli- 
gible effect on the mass, which neither improves nor dimin- 
ishes the level of agreement between rescaled and reference 
masses, independently of the definition of concentration ap- 
plied. 

Regarding the recovery of the concentrations of the 
haloes in the target cosmology, the individual concentra- 
tions defined affected by a bias of ~ 10 
percent with a strong dependence with the concentration 
in the reference cosmology. On the other hand, the concen- 
trations computed u sing the mass dependence reported by 
iBullock et all (|200ir i show a negligible difference between 
rescaled and reference co smologies, with no dependence on 
the halo virial mass. The lZhao et al.l (|2008l ) concentrations 
are obtained using their dependence on the cosmological 
time when the main progenitor acquired 4 percent of the 
final FoF halo mass. This quantity can be computed using 
the information in the halo merger trees; however, due to 
the numerical resolution of our simulation, this quantity can 
only be calculated for the ~ 2000 most massive haloes. For 
these haloes, the concentration shows no offset between the 
rescaled and target simulat i ons b ut its dispersion is larger 
than for the iBullock et al.l (1200 ll ) concentrations, particu- 
larly for low c values. 

Even though the recovery of individual properties shows 
some biases and non-negligible scatter, the quantification 
of these effects can be used to gauge what studies can be 
performed with rescaled simulations and the statistics that 
are expected to be affected by these uncertainties. 

Our findings on the accuracy of the properties of haloes 
are consistent with those found by AWIO, which indicates 
that applying the full method to particles or haloes produces 
only small differences. We confirm this estimate using our 
simulations, where we only find negligible differences in the 
resulting accuracy when using individual particles. 

3.2 Examples of halo evolution 

For a more detailed comparison between the rescaled and 
reference halo catalogues, we also compare the detailed his- 
tories of individual friends-of-friends (FoF) haloes. 
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Figure 4. Mass function for rescaled and reference FoF halo cat- 
alogues (solid and dashed lines, respectively) for different red- 
shifts (different colours as indicated in the figure key) . The lower 
sub-panel shows the relative difference between the rescaled and 
reference catalogues, in the same colour scheme. 
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Figure 5. Distribution of halo spin parameters for the rescaled 
(dashed lines) and reference (solid) catalogues. The mean value 
and standard deviation of the log-normal distributions are shown 
in the legend. 



We choose the first, lOth, lOOth and 1000th most mas- 
sive haloes and show them in Figure O The panels show 
from top to bottom, the accretion history, the accretion via 
mergers and via smooth, individual particle infall, the evo- 
lution of the dimensionless spin parameter, and the frac- 
tion of mass in substructures. The number of particles per 
halo decreases to the right, and range from FoF masses of 
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4.2x lO^^'h^^Mo to l.lx lO^^h^^M©. Regardless of the halo 
mass, the individual accretion histories, spin parameter evo- 
lution, and fraction of mass in substructures, are remarkably 
well recovered. The only noticeable differences are shown in 
the mass accretion via mergers and individual infall of parti- 
cles, which cancel out as when one is overestimated the other 
compensates, showing the effects from confusion between ac- 
cretion of small mass haloes and infalling field particles. This 
effect can be appreciated more clearly in the mass accretion 
history of the 10th largest halo in Figure [S] 

These results indicate that individual merger trees 
should be reasonably suitable for rescaling and later use for 
semi-analytic type galaxy formation modeling. The history 
of mass accretion via mergers would ensure a 10 percent ac- 
curacy in the population of satellite galaxies, whose added 
stellar masses should be even more precise (according to the 
well recovered fraction of mass in substructures), and the 
accurate evolution of the spin parameter should ensure rea- 
sonable estimates of galaxy disc sizes. 

3.3 Statistical properties of haloes and their 
growth histories 

If the rescaled simulation is used to produce statistics of 
large populations, either of dark-matter haloes or simulated 
galaxies obtained via semi-analytic techniques, it is neces- 
sary to estimate how they are affected by the rescaling pro- 
cess. 

We first study the mass function of dark-matter haloes. 
Figure U compares the mass functions of the rescaled and 
reference haloes (dashed and solid lines, respectively) for 
different redshifts from z = to z — 2. The agreement 
is excellent with only minor differences. The relative dif- 
ferences shown in the lower panel indicate that there is a 
mild, although clear trend that starts as an underestima- 
tion of abundances at z = but then tends to overestimate 
them by larger amounts at higher redshifts. At low redshift 
the rescaled catalogue shows a slight underestimation of the 
number density of haloes which is very mild at low masses 
M ~ 10^''h~^MQ at the 1 percent level, but increases to 
a 5 percent lower abundances at M ~ lO^'^'^h'^M©. The 
highest redshift shown in the figure indicates a fiat overesti- 
mation of abundances of a ~ 3 — 5 percent. 

These effects are small in comparison to the current pre- 
cision of measurements of the mass function in clusters of 
galaxies (|Gladders et al.lboOTi ') and ga laxy luminosity func - 
tions (for example in the SDSS, see iBlanton et al.l |2003| ) 
which, in simulations, are highly infiuenced by the under- 
lying mass funct ion of dark-matter haloes (see for instance 
ICole et aLllioOol ). However, the determination of a power 
spectrum f rom a mass function requires high precision mea- 
surements l|Sanchez. Padilla fc Lambasll200 ^ and the accu- 
racy of the rescaling technique could be important in this 
case. 

The spin parameter influences the re sulting properties 
of disc galaxies in semi-analytic models JCole et al.l I2OO0I : 
iLagos. Padilla fc Corall2009l : iTecce et~al]|2010l ). Therefore it 
is important to check whether the rescaling of a numerical 
simulation could produce biased distributions of spin param- 
eters. Figure [S] shows these distributions for the rescaled and 
reference catalogues (dashed and solid lines, respectively), 
and as can be seen, the distributions are consistent with one 



another. The distributions of disc sizes would therefore be 
expected to be reliable in the rescaled simulation. This is 
also in agreement with the evolution of the spins of individ- 
ual matched rescaled and reference haloes shown in Figure 
[S] which are consistent with one another at all times. 

Given that the stellar population of a galaxy is highly 
dependent on the mass accretion history of a halo, either via 
mergers or smooth infall, a statistical study of the merger 
trees of haloes can also be used to test the adequacy of the 
rescaling. Figure |6] shows the differences between rescaled 
and reference models for four particular characteristics of 
merger trees. The quantity shown on the top-left panel cor- 
responds to the time when half the total mass of a halo was 
formed in the tree t*" (i.e. including all the mass in satel- 
lites that will later merge with the central halo); by sub- 
tracting any cooling and star formation timescales involved, 
this would correspond to a stellar age for the final galaxy. 
As can be seen, the differences show a clear peak at OGyr 
and a median at ~ 0.19 Gyr, with a mean witdh of 0.84Gyr. 
On a 2: = galaxy population the former difference would 
produce little effects on the resulting galaxies, but the latter 
may spuriously broaden distributions of galaxy properties 
such as colours. 

The second statistics, the time since the halo in the 
main branch of the merger tree attained half of its final 
mass, (top-right panel of Figure |6}, would be related to 
the time of assembly of the mass of the final galaxy and 
shows an offset of ~ 0.19Gyr and a mean witdh of 0.86Gyr. 
Similar results are obtained for the time since the last major 
merger t""" (bottom-left) which should correlate with the 
time since the last starburst (the actual starburst should be 
more recent than this quantity when taking into account the 
dynamical friction that affect galaxies) , which shows a simi- 
lar offset of ~ 0.12Gyr on average, and a similar distribution 
width of 0.69Gyr. 

Since the timescales of last starburst, star formation, 
and stellar mass assembly are affected in a similar way, we 
do not expect important changes on the spectral shape of a 
galaxy obtained using rescaled merger-trees, with the excep- 
tion of the effect from a general shift in age towards older 
populations by ~ 0.2Gyr. The lower-right panel of Figure [S] 
shows the histogram of the difference between the number 
of major mergers undergone by the rescaled and reference 
haloes between 2 = 4 and 0. As can be seen, the distribution 
is symmetric with a clear peak at zero (which contains more 
than 50 percent of the sample), and a maximum difference of 
2 or more major mergers for about 10 percent of the haloes. 

AWIO presented comparisons between semi-analitic 
galaxies obtained from rescaled and direct simulations, 
showing an effect that may be the product of these offsets, 
the slightly fainter K-band magnitudes (by 0.1 mag) in their 
rescaled simulations. 

The study of the effect of the rescaling of merger 
trees on the inferred positions of haloes performed in 
Section [3] showed that these are not badly affected by 
the process. However, the correlation function could show 
changes due to the dispersion in the positions shown in 
the top-left panel of Figure [2] Figure [7] shows the resulting 
cross-correlation functions from the rescaled and reference 
haloes. For this test we choose cross-correlations over auto- 
correlation functions in order to increase t he signal-to-noise 
ratio of our measurements (see for instance lBornancini et aU 
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Figure 6. Distributions of differences between merger-tree properties for dark-matter haloes in the rescaled and reference modeh Top-left: 
distribution of the error in the time since the mass added over all the branches of the merger tree reached half the final mass. Top-right: 
same as top- left for the time since the halo in the main branch of the merger tree attained half its final mass. Bottom-left: same as 
previous panels for the time since the last major merger (with at least a 0.3 mass ratio). Bottom-right: differences between the number 
of major mergers for haloes which have undergone at least one major merger. The median m and mean width id of the distributions are 
given in the legend of each panel; we also show the average quantities in the reference cosmology. For the upper panels and the lower-left 
panel, the units are Gyr. 



I2OO6I : iLacerna fc Padiliall201ll '). Halos of different lower lim- 
its on FoF mass are used as centres, whereas all haloes 
in the simulations are used as tracers for these measure- 
ments. We choose three different lower limits on the mass of 
centre haloes {Mc > 10^^h"^MQ, Mc > lO"h"^M0 and 
Mc > 3 X 10^^h~^MQ) such that they bracket the non-linear 
mass for the cosmology of simulation B, M^i ~ lO^"^h~^M0, 
around which the bias factor shows a clear increase (for 
haloes with mass Mni , the bias factor is 6 = 1, for more 
details see for instance ISheth. Mo fc Tormen|[200ll '). 

As can be seen, at scales r > 0.5, 2.5, and 4h~^Mpc 
for masses Mc > 1 x 10^^ 1 x lO" and 3 x lO^h-^M©, 
respectively, the precision of the correlation function of the 
rescaled haloes is better than 5 percent (see the ratios on the 
lower sub-panel). Given the offset in the number of particles 
between the reference and rescaled haloes, we also tested 
whether using the same equivalent lower limit in number 
of particles improved this comparison, but find very similar 
results. This is also the case for the effect of the dispersion 
between rescaled and reference halo masses (number of par- 
ticles), since the resulting change is of only a few percent. 



which influences the clustering amplitude by factors below 
the offsets originating from the rescaling procedure. 

This level of precision is of the order of that obtained 
for lar ge surveys at z = such as the SDSS (see for in- 
stance IZehavi et al. 2004 ). Furthermore, other effe cts such 
as the assembly bias (Gao. Springel fc Whitell2005i ) are ex- 
pected to produce variations on the clustering of haloes 
and, consequently, galaxies to a 10 percent, which indicates 
that rescaled haloes could be used to test the detectabil- 
ity of this particular effect. However. IWu. Rozo fc WechsleJ 
(I2008i ) point out that the precisio n of future surveys such 



as the Dark Energy Survey (DES, iTucker et alll2007l) and 
the Large Synoptic Survey Telescope (LSST. lAbell et al.l 
2009) will be much better than that in the SDSS, and could 
irequire new runs with the desired cosmology rather than 
rescaled halo catalogues. 

As a final test, we study the changes introduced by the 
rescaling of haloes in the fraction of FoF particles in sub- 
structures (i.e., not restricting substructures to the sphere 
contained in one virial radius). Figure [8] shows this frac- 
tion as a function of the number of particles per halo, for 
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Figure 7. Cross correlation functions between haloes in the 
rescaled A (dashed lines) and reference B (solid lines) catalogues 
with different lower limits on FoF mass (different colours shown 
in the key) and the full population of haloes (with at least 10 
particles). The lower sub-panel shows the relative differences be- 
tween the correlation functions obtained from the rescaled and 
reference catalogues; the shaded regions show the uncertainties 
computed using the boostrap resampling technique. 

the rescaled and reference halo catalogues (lines and dots in 
blue and red, respectively). Both the median relation, and 
the envelope of 80 percent of the haloes show an excellent 
agreement between the two catalogues, with differences that 
only amount to a 10 percent. This statistics is related to 
different properties of the galaxies inhabiting these haloes. 
On the one hand, it would influence the variation of the 
amplitude of clustering as a function of subhalo mass (or 
galaxy luminosity), since the average host halo mass of a 
selection of subhaloes can be affected if this fraction is not 
accurately recovered. On the other hand, this result is in 
agreement with our earlier claim that we would not expect 
to find correlated changes in the spectra of semi-analytic 
galaxies obtained from the rescaled trees. This is due to the 
fact that the changes in the characteristic timescales of the 
merger trees are slightly biased but these biases go all in the 
same direction, otherwise we would expect to find shifts in 
the fraction of mass in subhaloes. This is also the cause for 
the fraction of mass in substructures to remain consistent; 
this can happen if the overall histories are consistent, even 
if slightly shifted. 



4 TESTS ON VARYING THE 

EXTRAPOLATION BASELINE 

The method of adjusting the cosmological parameters of a 
set of haloes and their assembly histories can be very use- 
ful for exploring a cosmological parameter space using halo 
properties from numerical simulations. In this section we 
investigate how the agreement between the rescaled and ref- 
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Figure 8. Dependence of the median (solid lines), 10 and 90 
percentiles (dashed lines) of the difference between the recovered 
and true fraction of mass in substructures as a function of the 
number of particles per halo. Blue and red lines correspond to the 
rescaled and reference halo catalogues, respectively. The uprising 
of fsub at small masses may correspond to a resolution effect. 



erence haloes degrades as the baseline of the extrapolation 
on the parameter space varies. 

We choose to vary only two cosmological parameters 
for this test, the matter density parameter Q.m, and the am- 
plitude of fluctuations erg. In this case we take as a start- 
ing point the parameters of simulation BZowQ which are 
in a greement with the latest constraints on the cosmol- 
ogy l|jarosik el al.|[201ol ). and explore the resulting rescaled 
haloes for the B0-2, Bo-i, Bo+i, B0+2 simulations, which 
maintain the same cosmological parameters as simulation 
Blow but vary Qrn between Qm — 2(j, and Qm + 2a with 
(J = 0.03, the simulations Bs-2, Bs-i, Bs+i, Bs+2, which 
with respect to Blow only change erg between the values 
(78 — 2a, and as + 2a with a = 0.03, and a simulation in 
which both, Q.m and erg are increased in one standard de- 
viation, Bo+is+i (details on the simulation parameters are 
shown in Tables 1 and 2). In all cases, the value of a is sim- 
ilar to the standard deviation in these parameters from the 
WMAP7 results. Since we are using simulations with small 
boxes for this test, we only apply the reduced version of the 
AWIO algorithm in this section. 

Figure[9]shows average and median quantities that sum- 
marise the accuracy of the rescaled haloes. We include the 
variation in the number of particles per halo averaged over 
all haloes with 200 to 1000 particles and the 10 and 90 
percentiles of the distribution (top-left), the average rela- 
tive variation in the cumulative mass function for masses 



^ The choice of small box and low resolution is justified by the 
precision of the rescaled haloes in the low resolution simulation 
Blow shown in Figure [2] and owes to the large number of simu- 
lations involved in this analysis. 
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in the range 11 < log-^Q{M /h~^ Mq) < 13 and the disper- 
sion of the distribution (top-right), the average variation in 
the cross-correlation function between haloes of FoF masses 
> lO^^h~^M0 against the full halo population, in the range 
of scales 1 < r /h~^Mpc< 10 and its dispersion (lower-left), 
and the median difference in the assembly ages of rescaled 
and reference haloes and the 10 and 90 percentiles of the 
distribution (lower-right). In all cases, the a;— axis shows the 
number of WMAP7 standard deviations from the Blow pa- 
rameters. As can be seen, the effect from changing the mat- 
ter density parameter is in all cases more important than 
varying the amplitude of fluctuations erg. We also show the 
resulting biases and uncertainties in the solid and open green 
symbols when changing both and erg by +la, which fol- 
lows the expected tendency and results in larger biases and 
dispersions than varying only one parameter at a time. 

These statistics are chosen since they are representative 
of important variations in the galaxy population that can 
be obtained via semi-analytic models (for instance via the 
abundance of haloes or the assembly time-scales), or vari- 
ations in constraints on cosmological models via clustering 
and abundance measurements. 

The biases in the number of particles per halo and the 
number of haloes shown in the top panels are correlated, and 
consist on underestimates (overestimates) for lower (higher) 
values of either Q,rn and erg. This is due to the fact that a 
bias in the number of particles, which does not depend on the 
halo mass as was seen in the comparison between simulations 
A and B, Alow and Blow and Abig and Bbig (cf. Figure [2l 
although an underestimation of the number of particles per 
halo is present), is reflected in the mass function. The change 
from simulation A to B is approximately equivalent to —Sa 
and +la in erg and Qrn, respectively, which would correspond 
to Vl'^ + 3^(7 ~ —3a in the top-left panel of Figure[9l where 
the number of particles per halo in rescaled catalogues is 
underestimated, showing the consistency of our analysis. 

The lower-left panel of Figure [9] shows that the ampli- 
tude of the correlation function shifts from an overestima- 
tion when the parameters are lowered, to an underestima- 
tion when going to higher values of either Qm or as (with 
a larger variation in the former). Finally, the asembly times 
shown on the lower-right show almost no biases, but a slight 
tendency towards underestimating this quantity for positive 
a values on either ag or Qrn (and stronger for the latter). 

The increase in number of particles per halo, the space 
density of haloes, and the decrease in the amplitude of the 
correlation function (and possibly in the assembly time- 
scales) responds to the decreasing values of the rescaling 
parameter s. This can be seen in Figure [10] where we show 
in a colour gradient that, for a fixed value of flm (for in- 
stance, that of the Blow simulation, shown by the vertical 
line), a higher final value of matter density parameter re- 
quires s < 1. This figure also serves the purpose of allowing 
to rescale the expected biases shown in Figure[5]to a different 
base cosmology. The white diagonal lines are contours of the 
equal s values that satisfy AQrn = 0.03 at the Blow cosmo- 
logical parameter set, and show that for lower values of base 
cosmology f2,„, the biases of Figure[9]would only be obtained 
for smaller values of Af2„i. As this result is independent of 
the base and desired values of erg, this shows that a grid of 
cosmological simulations would need to more densely cover 
low Sim values, roughly following AQm ~ 0.03 x (Sim /0. 27), 
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Figure 10. Rescaling parameter s (color gradient) as a func- 
tion of the starting cosmology parameter Qm (x-axis) and AQrn, 
needed to reach the desired cosmology (y-axis). The remaining 
cosmological parameters are those of the B simulation, with the 
exception of erg which is set to 0.82 for the starting cosmology 
and 0.73 for the desired one. The white vertical line indicates 
the position of Qrn = 0.27, and the diagonal lines represent the 
contours crossing Q^n = 0.27 when AO^ = ±0.03. The scaling 
length parameter s is not sensitive to changes in erg. 

to ensure a stable accuracy in rescalings done along the Qrn 
parameter, for any value of erg. As for changes in erg, as 
this will entail changing the output redshift adopted for the 
rescaling rather than introducing a new dimension in the 
grid of simulations, this will not affect the design of the 
grid. 

These results can be used to set the maximum difference 
on either flm or erg that will be allowed when adopting the 
AWIO method to rescale dark-matter haloes from a given 
cosmology into different ones. This result is also expected to 
hold for larger simulation boxes as long as the full rescaling 
algorithm is applied. 



5 CONCLUSIONS 

We have tested the use of the full and a reduced form of the 
Angulo & White (2010) method to change the cosmology 
of a catalogue of simulated haloes with the aim to test its 
application to a cosmological parameter search using semi- 
analytic galaxies. Our tests on the effect of applying it to the 
haloes in a numerical simulation instead of to the individ- 
ual particles, shows a dramatic reduction of computational 
time. The main reason for this particular test is that future 
planned numerical cosmological simulations will be so large 
that the storage of their individual particles will be imprac- 
tical, in which case only halo properties will be available. 
Even in the event that the approach presented here is ap- 
plied to a simulation with available individual particles, the 
CPU time required to rescale DM haloes is orders of mag- 
nitude (at least one) smaller than it is required to rescale 
the particles. If one adds the time consumed in identifying 
DM haloes and constructing merger trees, the speed up of 
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Figure 9. Recovery of halo properties as the basehne for the extrapolation in cosmological parameters extends away from the base 
values. In all the panels the blue lines correspond to varying the erg parameter, and the red lines to variations in Vim- Solid lines show 
the mean (top-right and bottom-left) or median values (top-left and bottom-right). The dashed lines show the 10 and 90 percentiles of 
the resulting distributions in the upper-left and lower-right, and the dispersion in the other two panels. The green solid symbol shows 
the recovery of halo properties for the model in which both, f2m and trg are increased by Icr, shown at + I'^a on the x-axis; the 
open symbols show the 10 and 90 percentiles, or the dispersion. Top-left: variation in the cumulative number of haloes with masses 
11 < \o^-^^q(M / Mq) < 13. Top-right: variation in the number of particles per dark-matter halo, for haloes with 200 to 1000 particles. 
Bottom-left: variation in the correlation functions, averaged on 1 < r/h~^Mpc< 10. Bottom-right: difference between the assembly 
time-scales of rescaled and reference haloes. 



the process can reach an improvement of two or three or- 
ders of magnitude (R. Angulo, private communication). The 
reduced form of the rescaling algorithm further reduces the 
computational time, but at the expense of being limited to 
small simulation boxes. For the cosmological parameters of 
simulations A and B (see Tabled, the simulation box side 
should be ~ 50h~^Mpc or smaller. In such boxes, the cor- 
rection for the quasi-linear modes does not further improve 
the accuracy of the rescaling. 



We measured the achieved precision of the rescaled 
haloes extracted from a given set of cosmological param- 
eters (such as those adopted in the Millennium simulation) 
by comparing them to those extracted from numerical sim- 
ulations with the desired final parameters (in a first in- 
stance corresponding to a Bolshoi simulation cosmology), 
constructed using the same initial conditions as the simula- 
tion from which the rescaled haloes are taken. We compared 
individual properties, such as the number of particles per 
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halo, their pecuUar velocities (and variation in the direction 
of movement), concentration parameters, the mass of the 
haloes, and their positions. In all cases the precision of the 
recovered properties is comparable to what is obtained from 
the full modification of the individual particles in the simu- 
lation (AWIO). The level of precision is better than lOOkpc 
for the halo positions, and of 5 percent of their peculiar ve- 
locities. Both the number of particles per halo, and the virial 
halo mass are slightly underestimated, by 5 percent. 

It should be noticed, though, that in order to obtain this 
accuracy for the virial mass, it is necessary to perform an 
additional correction which takes into account the different 
values of Svir{z, rim) iu the rescaled and targe t cosmologies. 
The halo concentrations derived using the iBullock et al.l 
l|200ll ) recipe shows the lowest biases between rescaled and 
reference cosmologies, with no dependence on halo mass. 

We tested the effects from lowering the resolution of the 
simulation by a factor of 8 in the number of particles (for 
the same total periodic volume), and from increasing by a 
factor of ~ 100 the volume of the simulation (for the same 
number of particles). We expect changes in the two cases due 
to different reasons. In the case of changing the resolution, 
all the other parameters of the rescaling of the haloes were 
held fixed. In the change of the simulated volume, the range 
of masses over which the rescaling function was minimised 
was shift one order of magnitude to higher masses (see Table 
[21 . Lowering the resolution or increasing the volume do not 
change the precision in the recovered number of particles as 
a function of the number of particles in the halo. This implies 
that a lower resolution results in a lower precision at a fixed 
halo mass. The lower resolution does not produce any other 
significant biases. The increase of the volume does not result 
in a higher uncertainty in the recovered positions since the 
quasi-linear correction of long wavelength contributions in 
the AWIO method properly takes into account the different 
quasi-linear theory modes in the two models. 

Furthermore, we also explored the differences arising in 
the detailed history of growth of haloes, and find only mild 
displacements in particular events such as major mergers. 
In order to do this we checked haloes corresponding to the 
1st., 10th, 100th and 1000th most massive objects in the 
simulations, and the resulting agreement is independent of 
the halo mass, at least down to masses M ~ 8x lO^^h~^M0. 
These results indicate that rescaled haloes and halo histories 
can even be used independently with little effects on the 
resulting galaxy population. 

A direct application of this method consists on chang- 
ing the cosmological background of semi-analytic galaxies. 
In order to assess the possible systematic biases that such 
a population of galaxies would present if it was based on 
a catalogue of rescaled haloes, we studied the offset in sev- 
eral measures of the halo merger histories, on the halo mass 
funcion and their clustering properties, on their spin pa- 
rameter distributions, and in the number of major mergers 
experienced by the haloes. Neither the mass function nor 
the correlation amplitudes show important variations (< 5 
percent) within the range of masses explored, which is also 
the case for the distribution of spin parameters or frequency 
of major mergers. The only possible important discrepancy 
comes from a systematic offset in the ages of merger trees 
(assembly, star formation and last starburst), which are all 
biased toward smaller ages. Given that these are correlated. 



the relative shapes of the spectral energy distributions of 
galaxies of different types will remain almost unchanged, 
but will result in a global change of the population toward 
slightly bluer colours. 

This method can effectively allow to sample the cos- 
mological parameter space using fully non-linear simula- 
tions. Even though semi-analytic models were designed to 
understand the processes driving galaxy formation and evo- 
lution, the possible dependences of galaxy properties on 
the cosmological parameters could also be used to im- 
pose constrains on the latter. To do this one possible ap- 
proach is to use Monte-Carlo Markov Chain or similar anal- 
yses of the cosmological para met er space (see for in stance 
iHarker. Cole fc Jenkin3l2007l and iBower etal]|2010D . This 
would produce chains of parameter sets on which, in prin- 
ciple, a new numerical simulation would need to be run, 
in which the detection of haloes and merger trees, and a 
semi-analytic model, would need to be performed and ap- 
plied. This process becomes much more efficient when using 
rescaled haloes, which avoid several time-consuming pro- 
cesses. The lowest computational cost is achieved by rescal- 
ing haloes in boxes small enough so as to allow the use of 
the reduced method. 

At this point it is important to note that even if there 
are small systematic effects on the scaled catalogues and 
merger trees, these are negligible in comparison to the uncer- 
tainties in the semi-analytic modeling. Furthermore, given 
the fiexibility of semi-analytic models, it is in principle pos- 
sible to diminish the difference between rescaled and direct 
simulations, but in a way that is not related to the scaling 
formalism; therefore, this particular ability does not help to 
improve the use of rescalings in a fast search for baryonic 
and cosmological parameters together. 

We showed that in terms of the uncertainties in the 
matter density parameter ilm and amplitude of fiuctuations 
parametrised by erg, the biases on different halo properties 
increase considerably more in the case of adopting more dif- 
ferent values of f2,„ than of erg. These constraints can be 
used to set limits on how far the extrapolation can be ex- 
tended in terms of these parameters; this way only a limited 
number of numerical simulations located in strategically se- 
lected grid points would need to be run in order to cover in 
a continuous way a wide cosmological parameter space. 
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